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Primitive path analyses of entanglements are performed over a wide range of chain lengths for both bead 
spring and atomistic polyethylene polymer melts. Estimators for the entanglement length A^'e which operate 
on results for a single chain length N are shown to produce systematic 0{1/N) errors. The mathematical 
roots of these errors are identified as (a) treating chain ends as entanglements and (b) neglecting non-Gaussian 
corrections to chain and primitive path dimensions. The prefactors for the 0{1/N) errors may be large; in 
general their magnitude depends both on the polymer model and the method used to obtain primitive paths. We 
propose, derive and test new estimators which eliminate these systematic errors using information obtainable 
from the variation of entanglement characteristics with chain length. The new estimators produce accurate 
results for Ne from marginally entangled systems. Formulas based on direct enumeration of entanglements 
appear to converge faster and are simpler to apply. 
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I. INTRODUCTION 
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The features of polymer melt rheology are determined pri- 
marily by the random-walk-like structure of the constituent 
chains and the fact that chains cannot cross. The motion of 
sufficiently long chains is limited by "entanglements" which 
are topological constraints imposed by the other chains. These 
become important and dramatically change many melt proper- 
ties (e. g. diffusivity and viscosity) as the degree of polymer- 
ization becomes larger than the "entanglement length" A^e- 
The value of is both a key quantity measured in mechan- 
ical and rheological experiments and a key parameter in tube 
theories of dense polymeric systems [ 1 ] . 

Ne is usually considered to be a number set by chemistry 
and thermodynamic conditions (e. g. chain stiffness, concen- 
tration, and temperature). It has been empirically related to a 
"packing" length Ne oc {pb^)~^ ||3D, where p is monomer 
number density and &^ = {R'^^f {N — 1)) is the statistical seg- 
ment length of chains with end-to-end distance i?ee and mean 
degree of polymerization N. In terms of individual entan- 
glements, Ne is defined as the ratio between N and the mean 
number of entanglements per chain {Z),in the limit of infinite 
chain length. 



N 

Ne = Van -TTTT- 

N^oo {Z) 



(1) 



We call a function We (N) an A'^e-estimate if it has the prop- 
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N 



lim We (TV) = iVe 



(2) 



where A^e is a system dependent but A^-independent quan- 
tity. Comparing Eq. ([1]) with (|2| does not imply choosing 
Afe{N) = N/{Z). The typical experimental A^e-estimate 
uses the plateau modulus G% ul]: 



We (AT) = 



(3) 



where m is monomer mass, fee is the Boltzmann constant, and 
T is temperature. 

A closely related theoretical construct is the primitive path 
(PP), defined by Edwards [4] as the shortest path a chain fixed 
at its ends can follow without crossing any other chains. Ru- 
binstein and Helfand |5] realized that the entanglement net- 
work of a system could be obtained by reducing all chains 
to their PPs simultaneously. Such a reduction process is an- 
alytically intractable, but has recently been achieved by com- 
puter simulations |6, 7, 8, 9, 10, lll[lJ,[i3[] which generate 
networks of PPs from model polymer melts, glasses, random 
jammed packings and solutions. These simulations estimate 
A^e either from the chain statistics of the PPs |6, 8, 9] or from 
direct enumeration of entanglements (contacts between PPs) 
, which determines (Z). 

Chain-statistical and direct enumeration approaches pro- 
duce different results for A^g for the same atomistic configura- 
tions, suggesting that "rheological" and "topological" entan- 
glements are not equivalent lUSll . This discrepancy has been 
attributed to the fact that chemical distances between entan- 
glements are not uniform, but rather are drawn from broad 
distributions , even at equilibrium. Studies 

of how entanglement properties change with A^ are therefore 
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of obvious interest. Moreover, primitive path statistics enter 
recently developed sliplink-based models I ItIi . 

In this paper we seek an "ideal" iVg-estimate which ap- 
proaches iVg at the smallest possible N. There have been 
several attempts in the literature, summarized below, to derive 
iVe-estimates, but these have all exhibited poor convergence 
(i. e. by approaching only at large N ^ N^). Molecu- 
lar dynamics simulation times increase with chain length N 
approximately as at large N (relaxation time r cx N^-^ 
times system size oc N^^^), so improved A^g-estimates have 
obvious benefits for computationally efficient determination 
of Ne- By analyzing a number of coarse-grained and atom- 
istic systems, we find a rather general solution to this problem 
of setting up a iVg-estimator which allows to predict from 
weakly entangled linear polymer melts. 

The organization of this paper is as follows. Section HI] 
presents the polymer models used here and the topological 
analysis methods which provide us with the entanglement 
network (primitive paths). Section Hill distinguishes between 
valid and quickly converging (ideal) A^g-estimators, and dis- 
cusses some model- and method-independent issues with ex- 
isting estimators. Examples are given which highlight system- 
atic errors caused by improper treatment of chain ends and of 
the non-Gaussian statistics of chains and primitive paths. Sec- 
tion |IV] derives two (potentially) near-ideal estimators which 
extract from the variation of entanglement characteristics 
with N. Section [V] presents and discusses numerical results 
for these estimators for two very different model polymers. 
We verify that they are basically ideal, explain why this is so, 
and derive simplified forms which further illustrate the con- 
nection of to chain structure and entanglement statistics 
and are also near-ideal. Section|VT]contains conclusions, and 
two Appendices provide additional technical details. 



n. POLYMER MODELS AND METHODS 
A. Polymer model systems 

We have created thoroughly equilibrated configurations for 
two very different (but commonly used) model polymer melts; 
monodisperse 'Kremer-Grest' bead-spring chains, and atom- 
istic, polydisperse polyethylene. These two are chosen be- 
cause they have similar values of but very different chain 
stiffness constants C(oo). Polyethylene is much more "tightly 
entangled" [18] in the sense of having a much lower value of 
iVe/C(oo); cf. Tabs.|II]and|in] 

The bead spring model lfl9ll captures the features of poly- 
mers which are key to entanglement physics, most impor- 
tantly chain connectivity/uncrossability. Each chain contains 
N beads of mass m. All beads interact via the truncated 
and shifted Lennard- Jones potential Uu{r) = 4eLj[(o'/r)^^ — 
(a/rf - (a/rcY^ + [cr/rcf], where = is the cutoff 

radius and Uu{r) = for r > Vc- Here a is the bead diam- 
eter and elj is the binding energy, which are both set to 1; all 
quantities will thus be dimensionless and given in the conven- 
tional Lennard-Jones units. Covalent bonds between adjacent 
monomers on a chain are modeled using the FENE potential 



U{r) = — iA:i?Q ln[l — (r/Ro)'^], with the canonical parame- 
ter choices i?o = 1.5 and A: = 30 [19]. The equilibrium bond 
length is Iq w 0.96. This model is hereafter referred to as the 
"LJ + FENE" model. 

Values of the density and temperature (p = 0.85 and 
T = 1.0) are those typically used for melt simulations 
lllllll9'l. All systems contain 280,000 total beads. While all 
are monodisperse, we employ a wide range of chain lengths, 
4 < iV < 3500. Those with TV > 100 are equilibrated 
using the "double-bridging hybrid" (DBH) algorithm |20]. 
DBH uses molecular dynamics to update monomer positions 
and Monte Carlo chain-topology-altering moves 12111 to over- 
come the slow diffusive dynamics 1 1] of entangled chains. All 
equilibration simulations were performed using the LAMMPS 
|22I] molecular dynamics code. Ref. [19] predicted A^e ~ 35 
at the above-mentioned state point using various "rheolog- 
ical" measures applied to systems with N < 400, while a 
similar analysis in [23] predicted Ne ~ 75. 

In all simulations of the atomistic polyethylene (PE) melt, 
the united atom (UA) representation is adopted. Accordingly, 
carbon atoms along with their bonded hydrogens are lumped 
into single spherical interacting sites. There is no distinction 
between methyl and methylene units in the interaction poten- 
tials. All bond lengths are kept constant (^q = 1.54 A), while 
bending and torsion angles are respectivelygoverned by har- 
monic and sum-of-cosines potentials I24l l25ll . Pair interac- 
tions between all intermolecular neighbors, and intramolec- 
ular neighbors separated by more than three bonds, are de- 
scribed by the 12-6 Lennard-Jones potential. The parameters 
of the mathematical formulas for the bonded and non-bonded 
interactions are given in Refs. ifTol l2ll l24l l25ll . These inter- 
action potentials yield accurate predictions of the volumetric, 
structural and conformational properties of PEjnelts over a 
wide range of chain lengths and temperatures i2l[|24ll . 

All atomistic PE systems were equilibrated through 
Monte Carlo (MC) simulations based on advanced chain- 
connectivity-altering algorithms: the end-bridging [26] and 
double bridging [21, 25] moves along with their intramolecu- 
lar variants. The simulated systems are characterized by av- 
erage chain lengths from = 24 up to = 1000, with a 
small degree of polydispersity. Chain lengths are uniformly 
distributed over the interval [(1 - A)Ar, (1 + A)Af]. Here A, 
the half width of the uniform chain length distribution reduced 
by A^, is 0.5 and 0.4 for 24 < Af < 224 and 270 < A^ < 1000, 
respectively. More details about the MC scheme, including a 
full list of moves, attempt probabilities and acceptance rates, 
can be found elsewhere [24]. Equilibration at all length scales, 
which is essential to obtaining meaningful results from en- 
tanglement analyses [27], was verified using several metrics 
II24I1 . In this study, results are presented for T — 400 K and 
T = 450 K, both for P = 1 atm. 



B. Entanglement network and primitive paths 

For the melt configurations the reduction to primitive paths 
was performed using two methods, PPA and Z, using the pro- 
cedures described in Refs. ll6l[7il9l [T4il . PPA simulations used 
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TABLE I: Chain and Primitive Patli dimensions for PPA and Zl as well as number of kinks (Z) for Zl for the LJ + FENE polymer melt. All 
quantities given in reduced LJ units. It is remarkable that values obtained via Zl and PPA are very comparable, suggesting that chain thickness 
and slippage effects seem to cancel as discussed in (3]. 



LAMMPS and Z simulations used the Zl code lUS]. Both 
PPA and Zl analyses are performed for the LJ+FENE model, 
while only Zl analysis is performed for PE. In both meth- 
ods, all chain ends are fixed in space. Intrachain excluded 
volume interactions are disabled while chain uncrossability 
is retained. Both classical PPA 16] and geometrical meth- 
ods (Zl H or CReTA |0, uM) provide the configura- 
tion of the entanglement network and the contour lengths Lpp 
of each primitive path. In PPA, disabling intrachain excluded 
volume produces a tensile force |29] in chains which reduces 
the contour lengths. In Zl, contour lengths are monotonically 
reduced through geometrical moves in the limit of zero prim- 
itive chain thickness. In addition to Lpp and the configuration 
of the entanglement network, Zl analysis also yields the num- 
ber of interior "kinks" in the three-dimensional prim- 
itive path of each chain. (Z) is considered to be proportional 
to the number of entanglements, regardless of the details of 
the definition used to define an entanglement. 

Runs end when the mean length of the primitive paths, 
(Lpp), and/or the mean number of interior kinks per chain, 
(Z), converge. Self entanglements are neglected, but their 
number is inconsequential for the systems considered here 

. The CReTA method works similarly, and the conclusions 
reached here for Zl analysis should apply similarly to CReTA 
results iillil. 

Table U summarizes chain and primitive path dimensions 
as well as (Z) for LJh-FENE chains with 20 < iV < 3500. 
Statistically independent initial states were used so that the 
random error on all quantities is < 2.5%. It is remarkable 
that PPA and Zl data for (Lpp) and also (Lpp) are so sim- 
ilar, considering the differences between the contour length 
reduction methods. Relative to Zl results, PPA values of 
(Lpp) are increased by finite chain thickness effects ifTll [30ll 



and decreased by chain end slipoff ifTill . Both these ef- 
fects should decrease in strength as N increases, and indeed 
(Lpp)ppA/(Lpp)zi decreases from ^ 1.17 to ~ 1.06 over the 
range 20 < iV < 3500. A very comparable trend is offered 

by (Lpp)2p^/(L 

pp/zi- 

PPA results for the shortest chains {N < 20) are not pre- 
sented. Standard PPA is unreliable for very short chains be- 
cause the presence of a high concentration of fixed chain ends 
combined with the finite bead diameter effectively inhibits re- 
laxation ifrn [30I1 . These problems are even worse for topo- 
logical analysis of lattice polymer systems, see e. g. Ref. (TjO- 
In the following, where (as will be shown) accurate data from 
very short chains is important, we focus on Zl results. 



III. TOWARDS VALID ESTIMATORS 

A basic task of topological analysis is to calculate from 
the full microscopic configuration of the entanglement net- 
work. The simplest approaches employ only the mean square 
end-to-end distance of chains, (i?ee) and either the mean 
length of the primitive paths, (Lpp) or the mean number of 
kinks, (Z). Notice that (Z) is not an integer, but semiposi- 
tive, (Z) > 0. In order to estimate from weakly entangled 
systems one of course needs physical insight; when this is 
limited, a good iVg-estimator can only be guessed. 

Some restrictions arise from a purely mathematical view- 
point. A valid estimator A/^(Af) has the following properties. 
It 

(i) obeys Eq. ^ and uses information from polymer con- 
figurations whose mean chain length does not exceed 

TV; 
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(ii) either yields Me (N) > N or leaves TVe {N) undefined 
for a system of completely unentangled {{Z) — 0) 
chains. 

An "ideal" estimator we define to 

(iii) correctly predict for all N exceeding N^, or for all 
(Z) exceeding unity. 

Accordingly, for an ideal estimator, the following weaker con- 
ditions hold. An ideal estimator 

(iv) diverges for a system of rodlike chains possessing = 
oo, and 

(v) exhibits Afe (N) < N when each chain has in average 
more than a single entanglement, (Z) > 1. 

The following two subsections repeat earlier approaches to 
estimate N^. Basic considerations of finite chain length ef- 
fects, errors from improper treatment of non-Gaussian struc- 
ture, and the general behavior of quantities entering are 
discussed. These subsections are meant to prepare the reader 
for the ideal estimators to be presented in Sec. |IV] They re- 
flect the chronology of our search for better estimators and 
help the reader to understand the magnitude of improvements 
presented in Section [V] The arguments given here ultimately 
point the way to construct ideal estimators. 



A. Non-Ideal Estimators 

Modelling primitive paths as random walks, Everaers et. al. 
10] developed an estimator (which we denote as "classical S- 
coil") which operates on results for configurations ("coils") of 
a single (S) chain length, 



Ue{N)^{N~l) 



{L 



(4) 



pp/ 



The classical S-coil estimate dU is useful because (for long 
chains) it relates changes in chain structure to rheological 
trends |0, [3. However, while it fulfills basic requirements 
(i) and (ii) (both unentangled and rodlike chains have i?ee = 
ipp), it lacks properties (iii) and (iv). As the exact relaction 
of (Ree) / {Lpp}'^ and (Z) is unknown, it is a priori unclear 
wether it has property (v). 

The corresponding estimator operating on the number of 
kinks, (Z), and originally employed in [7'], denoted here as 
"classical S-kink", is 



N{N - 1) 
{Z){N -1) + N' 



(5) 



which fulfills the basic requirements (i) and (ii), and also (v), 
but lacks (iii) according to Ref. [ 10] and (iv) by definition. The 
presence of both A^ — 1 and N in Eqs. (|5]), and subsequent 
estimators reflects the fact that it is the existence of a bond 
rather than a bead which is responsible for the presence or 
absence of an entanglement between two chain contours. 
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FIG. 1: (color online) Performance of classical A'^e -estimators 
Me{N) based on coils Q (upper two curves) and kinks l|5j. Data is 
shown for the two model polymer melts studied in this manuscript. 
The trends with A'^ is in agreement with published results for other 
systems The convergence behavior is poor, 

as Ne = limN ^oa J^e{N) obviously cannot be extrapolated study- 
ing chains with A'^ < 100, while A^e turns out to stay well below 100 
for both systems. An "ideal estimator", as defined in Sec. [nil would 
converge when A' exceeds Ne or earlier. 



The performance of the two classical estimators ^ and 
(|5]l for the two polymer models considered here is illustrated 
in Fig. [T] Values of Afe{N) converge very slowly with in- 
creasing N. As expected from their form, but contrary to 
both rheological intuition and condition (ii) , values of Af^ (N) 
drop strongly with decreasing N. For marginally entangled 
chains (where N is just large enough so that (Z) is small but 
nonzero), both classical estimators yield A^e (iV) < A^ — 1. For 
example, for N = 20, they both predict We (iV) = 17, which 
is close to the (improper) upper bound A^ — 1 = 19. This pre- 
diction obviously has no connection to the actual topology of 
the system. 

Thus Eqs. ^ and Q always underestimate, but never over- 
estimate Ne- This feature of the two estimators in the limit of 
unentangled chains is particularly (if retrospectively) disap- 
pointing, as it is incompatible with goal (iii). Similar behav- 
ior was reported (but not analyzed as in this paper) in Refs. 

Other previously published iVg-estimators ^ [Til, [3 l34ll 
also have some, but not all, of properties (i)-(v). One of the 
most promising was proposed in Ref. i^. It estimates A^e 
from the internal statistics of primitive paths, for a single A^. 
The squared Euclidean distances {R^{n)) between monomers 
separated by chemical distance n < A^ — 1 after topological 
analysis (i. e., the chain statistics of the primitive paths) were 
fit [9] to those of a freely rotating chain with fixed bond length 
fixed bending angle. was then identified with the chain 
stiffness constant C(oo) of the freely rotating chain 135]. This 
estimator does not obviously fail to meet any of conditions 
(i)-(v). In Ref. ]9] it gave values of Afe{N) which decreased 
more slowly than Eq. ^ as N decreased. Unfortunately, its 
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predictions agree with Eq. (|4|l at moderate N > 100 and thus 
it fails condition (iii). 

New S-estimators based on modifications to Eqs. (01) and ^ 
may be proposed. During the course of developing ideal esti- 
mators (to be introduced in Sec. llVb . we developed two mod- 
ified single chain length estimators which tend to approach 
Ne from above rather from below. These are the "modified 
S-kink" estimator 



160r 



N 



(6) 



and the mathematically similar "modified S-coil" estimator 

\ -1 



K{N) = {N-1) 



1 



(7) 



A motivation for the use of (L^) rather than (ipp)^ in Eq. (|7]i 
appears in Appendix lAl Fig.|2lshows results for Eqs. ^ and 
(|7]l for the same systems analyzed in Fig. [U Both modified 
single-chain estimators give J\fe{N) = oo for unentangled 
chains, thus fulfilling criterion (iv) in addition to (i) and (ii), 
but they still fail to fulfill goal (iii) since they tend to overesti- 
mate for weakly entangled chains. 



Errors from improper treatment of non-Gaussian structure 
and cliain ends 



Critically, none of the above-mentioned estimators seem to 
be able to predict Afe (N) — Ne for weakly entangled systems 
with a slightly positive (Z) < 1. All above-cited previous 
works as well as Eq. ^ have only produced convergence for 
N ^ Ne, and we are not aware of any studies where conver- 
gence has been achieved at « Ne, i. e., we are not aware 
of the former existence of any ideal A'^e-estimator. However, 
the failure of so many previous attempts both makes it worth 
examining the common reasons why they have failed, and in 
fact points the way to creating ideal iVg-estimators. 

To leading order in e = {N — 1)^^ (i. e. the inverse number 
of bonds), data for a wide variety of model polymers (see e.g. 
Refs. SIM HI), as well as the data obtained in this study 
(see Fig. [3^) are consistent with 



{Ri){e) = D/e-Y, 



(8) 



where the relative magnitudes of the constant coefficients Y 
and D depend on factors such as chain stiffness, molecular 
details and thermodynamic conditions. 

Also, orientations of successive PP segments are correlated 
mi], so (Lpp)^ should not be simply quadratic in chain length. 
The expected leading order behavior of (ipp)^ is 



(Lpp)2(e)=A/62 + S/6, 



(9) 



where B contains contributions from non-Gaussian statistics 
and contour length fluctuations 1 1]. Relationships (jS) and ^ 
are consistent with data reported elsewhere (e. g. Refs. 1241 
[32ll ) as well with our own data, as shown in Fig. 13 
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FIG. 2: (color online) (a) Performance of modified S-kink estima- 
tor ^ (lower two curves) and the modified S-coil Q (upper two 
curves), which approach from above. Data are for the same sys- 
tems analyzed in Fig. [3 The single-configuration estimator for kinks 
exhibits an improved convergence behavior compared with Un- 
der circumstances discussed in Sec.jV] application of both modified 
and original classical estimators allows one to obtain lower and up- 
per bounds on A'e which tighten with increasing TV. (b) Shown are 
the relative differences ("gap [%]") between A/'e( A'') values shown in 
Fig. [T] and the ones plotted in part (a) of the current graph. Differ- 
ences are smaller for A'^e estimated from kinks (lower two curves). 



At this point it is worthwhile to mention that we are going 
to make use of (|9]l, which is able to capture our results for 
(Lpp) down to chain lengths N small compared with Ne, to 
devise an ideal estimator in Sec.|lVl Relationship ^ however, 
as we will see, will not be required to hold to devise an ideal 
estimator 

Inserting Eqs. (|8]l and (|9]l into the classical and modified S- 
coil Eqs. (|4]l and O respectively give, to leading order in e. 



^,(„)i^-4I±££, + o(.=,, (10.) 
UAN)9 !^ + £l^^L^. + oi.% ,10b) 

Thus non-Gaussian structure of both chains and primitive 
paths naturally lead to systematic 0(e) ~ 0{1/N) errors in 
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ever, the source is chains being too short to be in the asymp- 
totic entangled limit defined by Eq. ([T). 

A key to understanding the failure of previous N^- 
estimators is that differences in the prefactors of the 0(e) 
errors [Eqs. (fTOl i and ( fT2] l1 arise from different treatment of 
chain ends. The classical S-kink equation ^ underestimates 
Ne as long as G + H < 1, and the modified S-kink equation 
(|6]l strictly overestimates N^, since both G and H are positive. 
Similarly, the prefactor (AY + BD- D'^)/A^ (Eq.[I^ of the 
systematic 0(e) error in the modified S-coil equation Q con- 
tains two contributions of different origins. [AY + BD) jA? 
arises from the Gaussian-chain approximation used, while 
—D'^/A^ arises from the attempt to correct for chain ends ef- 
fects (i. e. the "-1"). 

We have determined the coefficients A, B, D, G, H, and 
Y using all available data from our simulations; their values 
for both polymer models are shown in Tab. Coinciden- 
tally, for LJh-FENE chains, [ AY + BD)/A^ ~ 5 x 10^ and 
D'^ /A^ ~ 7 X 10"^. The systematic 0(e) eiTor for the mod- 
ified S-coil (|7]l is actually small for LJh-FENE systems due 
to the near-cancellation of its contributing terms. There is no 
reason to believe this behavior is general, and tests on addi- 
tional polymer models would be necessary i36tl to better char- 
acterize how rapidly the modified S-coil typically converges. 
However, it is reasonable to expect it typically converges more 
rapidly than the classical S-coil (|4|. 

Before turning to ideal estimators, we mention that the 
modified S-kink ^ can be regarded as corrected version of 
classical S-kink Q, as it eliminates an 0(e) error from the 
latter, and thus converges faster. 



FIG. 3: (a) Testing the applicability of Eq. ((S} which predicts linear 
behavior (slope oc D, offset oc Y) in this representation. We obtain 
C(oo) ^ 1.85 and C(oo) ^ 8.3 for the LJ+FENE and PE models, 
respectively (cf. Tab. Hit, (b) Testing the validity of Eq. ^ for both 
types of melts (slope oc A, offset oc B). The linear relationship is 
employed to derive estimator l ll5t in Sec. lIVBl (a-l-b) Data for larger 
TV are not shown but also agree to all displayed fit lines, to within 
statistical errors. 

Similarly, {Z) necessarily scales as e^^ in the N ^ oo 
limit. In the same spirit as the above analysis, and noting the 
failure (Fig.|2]i of Eqs. (|5]l and (|6]l to meet condition (iii), let us 
hypothesize that finite chain length leads to the leading order 
behavior 



{Z){e)^G/e-H, 



(11) 



where G and H are both positive. This assumption is actually 
consistent with the data in Tab.Uand previous works ll24ll : see 
also Sec.|IV] The classical and modified S-kink Eqs. (|5]l and 
(|6]l then become 



illljnEl 1 

G 



1-G~H 
G2 



e + 0(e2), (12a) 



Again, systematic 0(e) errors are predicted. In this case, how- 



IV. IDEAL ESTIMATORS 

Given the prevalence of subtle systematic 0(e) errors in 
non-ideal A^g-estimators, it is reasonable to suppose that in 
developing an ideal estimator, one has the freedom to intro- 
duce system-dependent (but A^-independent) coefficients, e. 
g., c, c', and Zq, in equations for a valid JVe{N) such as 
N/{Z) + ce, N{1 - c'e)/(Z), or N/{{Z) + Zq). These 
formulae are all potentially valid estimators because they ful- 
fill the basic requirement (Eq. The coefficients are some- 
what related to each other, but have slightly different physical 
meanings. They fulfill conditions (i) and (ii) for arbitrary c 
and c', but only if Z^ < 1. Note that finding an ideal A'e- 
estimator neither depends on the interpretation of (Z) or re- 
quires a priori knowledge of the numerical values of the coef- 
ficients. However, these numerical values are required to turn 
the above three expressions into A'^g-estimates before they can 
be applied. As these numerical values are certainly sensitive 
to system features like chain thickness and stiffness, it is im- 
possible to determine them from a single set of (ipp), {Rgg), 
and (Z) values. 

The best possible estimator gives J\fe{N) = for (Z) ^ 
1, but such an estimator would have to rely on incomplete 
information, some model assumptions, or make use of some 
'universal' features of entangled systems such as those sug- 
gested by Refs. {l^, [Till . We make use of two such findings 
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lo 


C(oo) 


D 


Y 


A 


B 


G 


H 


system 






cf. Eq. ^ 


cf. Eq. 


cf. Eq, 


. o 


LJ + FENE 


0.964 


1.852 


1.72 


3.55 


0.020 


1.04 


0.020 


0.12 


polyethylene 450 K 


1.54 A 


8.318 


19.7 A^ 


131.4 A^ 


0.22 A^ 


8.58 A^ 


0.023 


0.20 


polyethylene 400 K 


1.54 A 


8.535 


20.2 A^ 


85.3 A^ 


0.24 A^ 


9.37 A^ 


0.025 


0.19 



TABLE II: Data obtained via Zl. The coefficients D, Y, A, and B have been obtained from a least square fit to the available data (covering 
A'^ 2> A^e) for (-Ree) and (Lpp), according to Eqs. ^ and Similarly, coefficients G and H derive from the measured {Z) via Jilt . 



(SecHnll: for the polymer models considered here, both (Z) 
and {Lp^y /{N — 1) are linear in N above certain characteris- 
tic thresholds. Further supporting data for atomistic polyethy- 
lene have been reported recently by some of us ^27] . 

For both models considered here, the 'characteristic thresh- 
olds' are located at (Z) < 1 and N < N^, allowing us to make 
use of the 'linearities' to construct ideal A^e-estimators. We 
now derive two near-ideal A'^e-estimators, for kinks and coils 
respectively. These estimators operate on multiple (M) sys- 
tems with different chain lengths, rather than on a single con- 
figuration, and will be denoted as M-coil and M-kink in order 
to clearly distinguish between S- and M-estimators. Careful 
empirical tests of the new estimators' validity is quite essen- 
tial, and will be given in Sec. |V] 

Below, the idea behind the different roles of Eqs. (O, (|9|l, 
and (fTTT i is that the statistics of the entanglement network 
can be expected to be decoupled from the fractal dimen- 
sion of the atomistic chain, because entanglements arise from 
inter-chain rather than intra-chain configurational properties. 
The estimator we develop in the following section will, in 
fact, potentially be applicable to non-Gaussian chains where 
(i?gg) cx e~'^ (with 1 < /i < 2), as well as less-flexible poly- 
mers (like actin [37] or dendronized polymers |38]) for which 
A^'e is 111 8.] of the order of a "persistence length" of the atom- 
istic chain. 



A. The M-kink estimator 



as "M-kink", of extraordinary simplicity: 



K{N) dN ' 

M-kink is strictly an ideal estimator (i. e. it satisfies all five 
conditions proposed in Section Ullb provided A^i < N^. It 
eliminates the unknown coefficient in the linear relationship, 
and identifies to be responsible for the ultimate slope of 
{Z){N). This is analogous with measurements of diffusion 
coefficients, where one eliminates ballistic and other contribu- 
tions by taking a derivative. Application of Eq. requires 
studying more than a single chain length, which renders our 
M-kink estimator qualitatively different from the S-kink esti- 
mators. Data for {Z){N) for both polymer models, shown in 
Fig- in demonstrates that {Z) in fact becomes linear in A^ for 
{Z) below unity IH HI, thus confirming iVi < A^e- This 
suggests that can be estimated using data for {Z) from 
chains of lengths even below A"e. 

The occurrence of a nonvanishing A^i is rooted in the fact 
that a minimum polymeric contour length (of the order of 2Tit 
with polymer thickness f , subsequently corrected by chemical 
details) is needed for geometrical reasons to form an entan- 
glement (or tight knot) |40]. This length {£) increases with the 
persistence length of the atomistic contour, and vanishes in the 
limit of infinitely thin polymers. This implies that determin- 
ing A'e from the slope we correct for a thickness effect, and 
A^^i is proportional to the thickness of the atomistic polymer. 



Beyond some a priori unknown chain length A^i, we know 
that {Z) (as determined via Zl or CReTA) varies linearly with 
N, i. e. {Z) = GN + Zo (with G > 0, and Zo = -{G+H) > 
— 1 in the notation of Eq. [TTJ. We recall that an ideal A^e- 
estimator implies, according to condition (iii), that 

(vi) dAfe{N) /dN = for A^ > A^i, and 

(vii) A^i < Ne 

are necessary to produce = J\fe{Ni). Uniquely, A^e = 
1/G and Afe{N) = A^e for all A^ > A^i. Using flie hn- 
ear relationship between (Z) and A^ we thus propose (a) 
Af^{N) = N/{{Z) - Zo), where Zo = Zq{N) is the coef- 
ficient determined from data collected up to chain length A^. 
Note that (a) is identical with the A^'g-estimator suggested on 
mathematical grounds at the beginning of this section. 

However, Afe{N) — 1/G can be equivalently obtained 
from (b) Afe{N) = dN/d{Z). This is an estimator, denoted 



B. The M-coil estimator 

Next, we motivate and derive a near-ideal estimator for 
use with coil-properties (i?^^) and (ipp) (obtained via PPA, 
CReTA or Zl). Flory's characteristic ratio C{N) is defined 
through the identity ll4lll42ll 

(Ri) ^{N~ l)llC{N). (14) 

Equation (fl4] i is exact by construction; the A^-dependence of 
C{N) characterizes the (non)-Gaussian structure of chains. 
In general, G{N) > 1 if A^ > 1. For (mathematically) ideal 
chains, G{N) is related to the persistence length Ip |35]. This 
allows the chain stiffness constant C(oo) = limjv^oo C{N) 
to be calculated from short chains for any sort of ideal chain, 
including random walks, freely rotating chains, wormlike 
chains, etc. Simulations on dense chain packings show [43l] 
that the value of C{oo) = 1.48 is a universal lower limit for 
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excluded volume, flexible chain molecules. For real chains 
like polyethylene, chains much longer than Ip need to be stud- 
ied to characterize C{N), cf. Ref. ll24ll . We assume knowl- 
edge of C{N) as function of N from the atomistic configura- 
tions. 

To proceed, we make use of our finding that (Lpp)^/(A^ — 
1) is linear in N above a certain characteristic No, before 
{Z){N) has reached unity, i.e., we assume < to de- |B] 
rive an ideal estimator (fTSl l. The linear relationship clearly 
holds for both polymer models considered here (Ref. [24], Ta- 
ble|l] Figs.[3]5 and|4]i, and has already been formulated in Eq. 

Next we relate (Lpp) and {Z) for large > iVe by 
a simple argument: the length of the primitive path, ipp, is 
ifisll the number of "entanglement nodes", N/N^, times the 
mean Euclidean distance £e between such nodes. This dis- 
tance (£e) equals the mean end-to-end distance of the atom- 
istic chain with monomers. We thus expect that up to 
a factor of order unity (related to fluctuations in le 03511 ). 



limjv^o,(Lpp)2 = {N/N,)HN, - l)llC{N,). 




FIG. 4: Zl results for the two model polymer melts. Testing the 
applicability of Eq. dl It which predicts linear behavior in this repre- 
sentation (slope G, offset H). Clearly (Z) (N) becomes linear at an 
N for which {Z} < 1. This implies iVi < and that Eq. l fT3ll can 
be an ideal estimator. An interpretation for A'^i is given in Sec. lIV Al 
Data for larger TV is not shown here, but the slope d{Z)/dN does 
not significantly change with increasing A^. 



By following the procedure of Section |IV A| we aiTive at 
an A^e-estimator, denoted as "M-coil", using coil properties 
alone: 



x=j\r^{N) 



d 

dN 



p2 



(15) 



where 



{N — 1)Iq, and C{x) is the characteristic ra- 
tio for a chain with J\fe{N) monomers. This estimator fulfills 
all conditions from our above definition of an ideal estimator 
As for M-kink, the derivative in the M-coil Eq. (flSl l signals 
that we have to measure (Lpp) as function of N rather than a 
single value to estimate N^- The convergence properties are 
not as clear a priori as they are for the M-kink estimator Eq. 



( HjI ). as this derivation required an approximation. In prac- 
tice, one must simulate systems with increasing N until the 
M-coil converges. There is no apparent way to come up with 
an TVe-estimator from coil quantities which converges before 
N reaches N^- This is a noticeable difference between the es- 
timators from coils and kinks (M-kink). Technical considera- 
tions in the application of Eq. ( fTSl l are discussed in Appendix 



V. NUMERICAL RESULTS AND DISCUSSION 



The data in Tab. U and a similar set for atomistic polyethy- 
lene (configurations from Ref. |24]), will now be used to test 
the M-estimators. Figure |5] shows results for the M-kink es- 
timator (Eq. [T3] l and M-coil estimator (Eq. [T5] l for the same 
systems analyzed in Figs. [T] and |2] Comparison of these fig- 
ures shows that the M- estimators indeed converge faster than 
the S- estimators (Eqs. |4}{7]). Moreover, comparison to Fig. 
ID shows that the M- estimators converge for marginally en- 
tangled systems; values of JVe{N) approach before (Z) 
far exceeds unity. These show that Eqs. ( fT3] l and ( fTSl ) are es- 
sentially "ideal", meeting all of conditions (i)-(v). The kink 
estimator performs slightly better, presumably because of the 
approximations made in deriving Eq. ( fTSl l. 

For < 50, values of J\fe{N) from M-kink ST3[ increase 
with decreasing A^. As shown in Fig. |5] Afe{N) appears to 
be diverging as A^ ^ 0. The precise nature of the divergence 
is unimportant. For example, A^ = 20 chains have (Z) = 
0.127, and the vast majority have zero entanglements, so the 
prediction A/'e (20) = 192 > 20 of modified S-kink ©just 
signals that we are deep in the unentangled regime, where A'e 
cannot yet be estimated. 

The fast convergence of the M-kink estimator can be better 
understood by plugging Eq. (fTTT i into M-kink ( fT3] ). This pro- 
duces a special case of the M-kink estimator, which is only 
asymptotically correct, and can be used when (fTTT l holds. We 
refer to it as the "approximate M-kink" estimator: 



1 

G' 



(16) 



Here, G is the coefficient in the linear relationship between 
{Z) and A^ obtained from data collected up to chain length A^, 
and thus Afe{N) depends on A^. Note that the derivative with 
respect to A^ in Eq. (fTSI ) removes the 0(e) errors! This is a 
major difference with respect to all S-estimators (the estima- 
tor used in |9] can be considered as intermediate between S- 
and M-estimators). 

In a similar attempt to rationalize the fast convergence of 
the M-coil estimator, we insert Eqs. ([8]l and ^ into Eq. ( fTSl ). 
This yields, accordingly, the "approximate M-coil" estimator. 



1 



D + VD^ - 4AY 
2A ■ 



(17) 



Like Eq. ( fTSI l. Eq. ( fTTb has no 0(e) corrections. Again, this 
arises from the "M" approach of taking derivatives with re- 
spect to A^. In both cases, the use of the derivatives re- 
moves undesirable effects related to proper treatment of chain 
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AAe(iVe) 




iVe A4(iVe) 








M-coil 


approximate M-coil 


simplified M-coil 


M-kink 


approximate M-kink 


system 


Eq. (dill 


Eq. Oil) 


Eq. (Ull 


Eq. (nil 


Eq. (dill 


LJ + FENE 


86.1 87.8 


85.1 89.6 


86.2 90.1 


48.9 46.3 


48.5 55.7 


polyethylene 450 K 


84.0 83.4 


84.4 84.5 


90.6 90.1 


44.2 42.2 


43.3 38.8 


polyethylene 400 K 


82.3 80.1 


80.5 77.8 


83.9 84.1 


41.5 38.5 


40.1 36.3 



TABLE III: Data obtained via Zl. Selected results for N^{N) for all near-ideal M-coil and M-kink estimators defined in this manuscript. 
For each estimator, two characteristic values are shown: TVe uses all available TV (up to A'^ = 3500 and A'^ = 1000 for the LJ+FENE and 
PE models, respectively), and Afe{Ne) uses only data from short chains with N < Ne (cf. Tab. U^. Values of Me{Ne) are thus obtained at 
moderate computational cost, and are all in overall agreement with N^.. Approximate M-coil (M-kink) results should coincide with M-coil 
(M-kink) results, if the relationships {Sj, ^ and dllt . respectively, accurately hold. The simplified M-coil does not take into account the effect 
of C{N). M-coil (M-kink) is the estimator with the least assumptions involved, if A'^e needs to be estimated from coil (kink) information 
(see also Appendix iBt. The fact that for all these estimators Afe{Ne) ~ A'^e gives sufficient evidence that these are in fact ideal estimators, in 
sharp contrast to most S-estimators, quantitatively discussed in Tab. lIVI Note that the very similar values of Ns reported for LJ+FENE and 
PE systems are a pure coincidence arising from their similar values of D/A (Tab.Ull cf. Eq.|18l) 





Ne Afe(N,) 


Ne AfeiNe) 


Ne Afe{Ne) 


Ne Afe{Ne) 




classical S-coil 


modified S-coil 


classical S-kink 


modified S-kink 


system 


Eq. iH 


Eq. 


Eq. g} 


Eq. 


LJ + FENE 


86.1 40.0 


86.1 129.7 


48.9 31.7 


48.9 51.2 


polyethylene 450 K 


84.0 39.7 


84.0 192.8 


44.2 30.4 


44.2 48.3 


polyethylene 400 K 


82.3 37.3 


82.3 191.8 


41.5 28.8 


41.5 44.9 



TABLE IV: Data obtained via Zl. For comparison with Tab. IIIII Performance of previous S-coil and S-kink estimators. Accurate A'e-values 
have been overtaken from M-coil and M-kink in Tab. IIIII Obviously, J\fe{Ne) is far from being close to A^e in all cases, while the deviations 
are strongest for the A'^e -estimates based on coils; the two kink measures seem to at least bracket the true A'^e (for the deeper reason that Zq, 
introduced in Sec. lIV Al must obey Zo £ [—1, 0]). 



ends. The approximate M-coil estimator is related only to 
the (in general, non-Gaussian) structure of chains and primi- 
tive paths. Finally, if the assumptions which lead to Eq. ( fTOl i 
hold, and in order to quantify the contributions to Eq. (ill , the 
above analysis combined with tube-theoretic considerations 
suggests another estimator, which we refer to as "simplified 
M-coil": 

^feiN) « ^. (18) 

The only dependence on N of the approximate and simpli- 
fied estimators, Eqs. (fT6Tl-(fT8]l, stems from the variation of A, 
D, G, and Y with N; these coefficients, which are obtained 
by linear interpolation, must generally be assumed be consid- 
ered to depend on the available range of studied chain lengths. 
When the variation of the coefficients is large, these three es- 
timators should not be used. 

Note that the simplified M-coil does not agree with M-coil 
if C{Ne) has not reached C{oo); though it may converge 
quickly, it cannot be ideal. For the systems under study, is 
large enough such that C{Ne) is quite close to C{oo) ||44|]. 
The simplified M-coil has a simple connection to polymer 
structure and the tube model fH). D — C(oo)/q = IqIk, 
where Ik is the Kuhn length |35]. The tube diameter (It is 
given by — I^IkN^, and hence A = (dT/Ne)'^. 

Table |III] quantifies the performance of the new M- 



estimators. The two presented values for each estimator 
J\fe{N) are the final N^, obtained by analyzing all available 
chain lengths, together with the value predicted by the estima- 
tor at = Ne (i. e. at the border between unentangled and 
entangled regimes, using only chains of length up to ~ Ne). 
For an ideal M-estimator these two numbers should be the 
same within statistical errors, here ~ 2.5%, and Ng should 
coincide with limjv^oo N/ (Z). All four M-estimators con- 
sidered here, the complete ones (Eqs. fTJlfTsT l as well as their 
approximate versions (Eq. [16] [T7] i satisfy these criteria. The 
simplified M-coil ( fTSl l is seen to converge quickly as well, 
but it does converge to an N^ which is above the one ob- 
tained via M-coil, because Y is positive (Y vanishes for an 
ideal random walk). TablejlVjshows corresponding results for 
the S-estimators, which all (as discussed above) are generally 
non-ideal. Still, the modified S-kink turns out to perform very 
well, simply because G + H 1 for our model systems, cf. 
Tab. In] 

For the LJ+FENE model, while the classical S-coil estima- 
tor (Eq.lU produces values of N^ consistent with published f% 
results, i. e. N^ ~ 70 for N = 350 and 500, values for these 
estimates based on the near-ideal M-estimators (cf. Tab. HIIl l 
and also the modified S-coil ^ rise above 80 for the longest 
chains considered here. The M-estimators based on chain and 
(Zl) primitive path dimensions converge to the value N^ ~ 85 
in the mostly unentangled regime, cf. Tab. |III] Thus all data 
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FIG. 5: (color online) (a) Performance of proposed estimators M- 
kink l ll3t (lower two curves with large symbols) and M-coil l llSt 
(upper two curves with large symbols); see also Appendix IbI Data 
is for the same systems analyzed in Figs. [T] and |2] Clearly, Afe{N) 
has converged for TV <C 100, and as shown by comparison to Fig. [4] 
J^e{N) approaches A'^e before (Z) exceeds unity. This allows us to 
estimate A'^^ from mostly unentangled systems, (b) Same data as in 
(a) vs. logio A'^, which allows the full range of A*' to be presented. For 
comparison, blue broken and red dashed lines for PE and LJ+FENE, 
respectively, show reference data for S-estimators, already presented 
in Figs.Eand|2] 



suggest that the "best" estimate of the entanglement length 
for flexible chains is well above the previously reported value. 
This is significant e. g, for quantifying the ratio N^/Nc, where 
Nc is the rheological crossover chain length where zero shear 
viscosity changes its scaling behavior from Rouse to reptation, 
and has been estimated as ~ 100 filH]. 

One could imagine fitting the squared contour length 
(Lpp(n)) of primitive path subsections ll46ll to {L^^{ri)) = 
An^ + Cn. and attempting to calculate Afe{N) = D/Ahy 
also fitting to {R^{n)) = Dn — Y , or developing other im- 
proved estimators for based on {R^^in)) and (ipp(n)). 
However, analysis along these lines failed to produce any es- 
timators better than those described above. In particular, no 
improvement over the method of Ref. ^ was found. 

It is important to notice that our Eq. dTsl l is not compati- 
ble with some earlier definitions of {Z) from coil quantities. 



because of the prefactor C {oo) / C {N e) ■ This prefactor had 
usually been omitted or not mentioned, since random walk 
statistics were clearly a convincing starting point. Assuming 
Gaussian statistics (constant C{N) for all N) hence under- 
estimates values of calculated from coil properties. This 
issue is also one of the reasons why the estimates between 
PPA and geometrical approaches differ. Another reason is 
given in [29]. Ratios between 1.3 and 2.5 between N,, calcu- 
lated from kinks and coils have been reported in the literature 
. The presented data exhibits ratios between 1.6 
and 2. A third reason that they differ is rooted in the fact that 
{Z) is not |7] uniquely defined from a given shortest, piece- 
wise straight path, as it is returned by Zl or CReTA. This 
additional discrepancy can only be resolved by matching re- 
sults for from kinks and coils, and by comparison with 
experiments. 

The classical S-kink (|5]l strictly underestimates and the 
modified S-kink ^ strictly overestimates (since both G 
and H are positive, and G -\- H < 1). 



VI. CONCLUSIONS 

Very significantly improved, near-ideal, and apparently 
polymer-model-independent estimators for were derived 
in this paper, M-coil (Eq. [15] to be used with PPA, Zl or 
CReTA) and M-kink (Eq. [13] Zl and CReTA only). They 
reduce, under further assumptions which seem valid for the 
model systems studied here, to approximate M-coil (Eq. [TTI l. 
simplified M-coil (Eq,[T8]l, and approximate M-kink (Eq.[T6]l. 
These estimators require simulation of multiple chain lengths, 
but have eliminated systematic 0(e) errors present in previous 
methods. This is important for the design of efficient simula- 
tion methods in the field of multiscale modeling of polymer 
melts. 

Furthermore, we have proposed variants of the original esti- 
mators. The two main problems with existing estimators were 
identified as: i) improper treatment of chain ends, and ii) non- 
treatment of the non-Gaussian statistics of chains and primi- 
tive paths ll36ll . Improper handling of thermal fluctuations was 
an additional problem relevant to very short chains. Issues i) 
and ii) lead to separate, independent 0(e) errors. Estimators 
based on direct enumeration of entanglements lack issue ii), 
and so are fundamentally advantageous for estimation of A^e- 
The new "M" estimators proposed here formally correct for 
the eiTors arising from effects i) and ii). The values of the M- 
coil and M-kink-estimators can be taken as "best estimates" 
for Ne when results are available for multiple chain lengths. 
The best estimator when only a single chain length is available 
is the modified S-kink, Eq. (O. 

We have shown that e(ipp)-^, (Z), and also {R"^^) are all 
linear in 1/e (thus linear in N) down to the mostly unentan- 
gled regime, and have used this information to derive the M- 
estimators and to improve the earlier ones. All coefficients 
in these linear relationships have been evaluated and listed 
in Tab. [n] The prefactors for the above mentioned 0(e) er- 
rors can be large, and depend both on the polymer model 
and method of topological analysis. These errors can produce 
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large changes in estimates of iVg for values of N typically 
considered in previous studies (e. g. Refs. ll6l [T2l[33ll ). This is 
significant in light of attempts to compare PPA results for Ng 
to values obtained by other methods |^ |23l [lH ll^l such as 
direct rheological measurement of the plateau modulus G%, 
evolution of the time-dependent structure factor S{q, t), and 
estimation of the disentanglement time cx 
Some conclusions of those studies may need to be reevaluated 
in light of the new data. 

The proposed M-estimators are to our knowledge the first 
estimators which exhibit all features required for an ideal esti- 
mator (a term which we made precise in Sec.|lll), and they have 
been physically motivated. They converge to iVg for weakly 
entangled systems (N < Nf.). They leave either undefined 
or infinite for rodlike chains (because C{N) = N for a rod). 
They predict Af^ (N) > N for a completely unentangled sys- 
tem, which is characterized by (Z) = and ipp — i?ee in 
accord with the definition of the primitive path which we have 
adopted in this work (see [29]). The appearance of the coef- 
ficient A^i suggests that there might be a minimum amount of 
material, A^i, needed to form a single entanglement (as ob- 
served for phantom chains |7]). If so, it can be expected to 
depend on the thickness of the atomistic chain and its stiff- 
ness as well as particle density. We expect our findings to be 
universal in the sense that they should apply to all sorts of real 
linear polymer chains in the melt state, and we have verified 
the assumptions underlying the M-estimators by direct com- 
parison with both atomistic semiflexible and coarse-grained 
flexible polymer melts. 

Refs. ifm [Tsll pointed out that primitive paths are not ran- 
dom walks, and that there appears to be more than one "topo- 
logical" entanglement per "rheological" entanglement; thus it 
is unsurprising that from coils is significantly larger than 
iVe from kinks (for details see Ref. 11811 ). The utility of any 
topological analysis of chains shorter than remains highly 
questionable, because the chains' dynamics are well described 
by the Rouse model iH, [l9t] and so they cannot be considered 
"fully entangled" in any meaningful way. However, it seems 
that the M-estimators developed in this work have the ability 
to extract information from a partial or even marginal degree 
of entanglement. 

The M-estimators could be applied in a post-processing 
step on existing configurations. For example, it should be of 
interest to study the effect of flow and deformation on entan- 
glement network characteristics in order to establish equations 
of motion for relevant coarse-grained variables characterizing 
the polymer melt. Shear and elongational flows have been 
studied for both polymer models considered here, but either 
Zl was not yet available at the time of these studies (4§\, or 
the chains were |49, 50] "too short", i. e. had (Z) <C 1. 

The apparent ability to accurately estimate even for 
weakly entangled systems may be useful for atomistic mod- 
els whose computational cost prohibits equilibrating large- 
systems, such as polymers containing bulky side groups. The 
procedure for removal of the 0{e) systematic errors, while 
clearly described here, requires performing analyses on a lim- 
ited number of configurations on a range of chain lengths, 
which is most easily undertaken for systems composed of 



"short, but not too short" chains. 
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APPENDIX A: TREATMENT OF THERMAL 
FLUCTUATIONS 

Ref. and other studies have typically used (Lpp) ^ rather 
than (ipp) in estimators for N^, such as the analogue for the 
modified S-coil O which reads 

AAe(A^) = (iV-l)(^^-l) \ (Al) 

However, Eq. dAlb gives pathological results for short chains 
due to existing thermal fluctuations of Lpp. Consider the un- 
entangled limit, where the entanglement density (denoted as 
Pe) vanishes. For an "ideal" topological analysis, ipp i?ee 
(from above) for each and every chain as pe ^ 0. How- 
ever, chain dimensions fluctuate in thermodynamic equilib- 
rium Ol. To leading order in the fluctuations, (Lpp)^ = 
(L^p) - (ALpp)2 = (Rl) - (Ai?ee)2, where A is "vari- 
ance of". So, even for an ideal topological analysis proce- 
dure, Eq. dAlb would predict a negative Me [N] — > — [N — 
l)(i?2j/(Ai?ee)^ as Pe ~^ 0. Negative (iV) are of course 
useless, but indeed, are predicted using our data in Tab. H] For 

= 20 (LJh-FENE melt), application of Eq. ( lAll ) yields neg- 
ative A/^(20). A term identical to the term in parenthesis in 
Eq. lAll was found to be negative for short chains in Ref. ll49ll . 
but was not used to directly calculate Afe{N) in their work, as 
its negative value was considered to signal (and to only occur 
in) the mostly unentangled regime. 

The reason to fix chain ends during PPA or Zl analysis is 
the assumption, implicit in Edwards' definition of the prim- 
itive path i^, that chains are entangled. In this context it is 
worthwhile mentioning that there are other definitions of PP's, 
for example one [51] where the length of the PP goes down to 
zero for the unentangled chain, and where chain ends are not 
fixed. 
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APPENDIX B: TECHNICAL CONSIDERATIONS IN USE OF 
THE M-COIL ESTIMATOR 
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FIG. 6: The graph demonstrates on how to graphically evaluate 
Afs{N) according to the M-coil estimator l ll5t . Shown are both the 
left hand side (Ihs) C{N) /N, and right hand side (rhs) of Eq. (15} for 
both types of polymer melts. The dotted red path makes an example 
on how to obtain J^^^N) ~ 87 for given — 48. Obviously, this 
value is quite identical with both, TV^ and J\fe(Ne), cf. Tab.lIIII The 
ratio C{N)/N (small points) monotonically decreases with increas- 
ing A'^, while the rhs (large symbols) reaches a plateau at the time 
A'^ has approached A^e (at the crosspoint), which is a distinguished 
feature of an ideal estimator. 

While the M-kink estimator (Eq.[T3]l is explicitly evaluated 
from the local derivative d{Z) /dN around N, our M-coil ex- 



pression, Eq. ( fTsT i, is only an implicit expression for the esti- 
mator Me (N). Formally, we need the inverse of C{N) /N to 
calculate Afe{N). In the following, we describe the procedure 
in order to prevent any ambiguities upon applying M-coil in 
practice. Fig. |6] shows both the left (Ihs) and right hand (rhs) 
sides of Eq. [13] versus N for our data. For any given N (say, 
= 48 for the PE data, where the dotted red line starts in 
Fig. |6]l, the 7V^(iV) estimate is the value at the ordinate for 
which the abscissa values for Ihs and rhs coincide (end of the 
red curve is at JVe (48) « 87). The same procedure is repeated 
for all N to arrive at Fig. |5] and particular values collected in 
the M-coil row of Tab. Hill The difference between Ihs and rhs 
can be used to estimate the difference between the largest N 
available and N^. If only short chains had been studied, only 
a part of this plot could have been drawn. 

Note that this procedure requires C{N) /N to be monoton- 
ically decreasing with N, and access to C{N) at sufficiently 
large N. While the former is essentially valid for all poly- 
mer models, the latter may pose a problem. Without reliable 
values for C{N) for N = Ng, there is no apparent way to 
come up with an M-coil which converges before reaches 
A^e- However, since C{N)/N decreases with increasing 
and ultimately reaches C{oo)/N behavior, in practice (and 
formally for ideal chains) C (N) can be estimated by extrapo- 
lation, and the necessary C{N) /N values could be added for 
chain lengths exceeding those studied. 

This issue disappears by construction when the largest sim- 
ulated A^ exceed JVe{N), so that the conditions for an ideal 
estimator are met in any case. Still, this is a noticeable and 
principal difference between the estimators from coils (M- 
coil) and kinks (M-kink). 
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